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An original method, developed by R. Kappus, aims at 
defining the inertial constants of an aircraft or a rocket, 
or any other structure, even flexible, without materializing 
any rotating axis. This paper presents a new synthesis of 
the method. The necessary equipment is very similar to that 
used normally for ground vibration tests. The only new de- 
vice to be provided is an elastic suspension for obtaining 
th a total natural modes corresponding to the movements of the 
structures as a solid. From the measurement of the general- 
ized masses of these modes it is possible to compute the in- 
ertial constants: center of inertia, tensor of inertia, mass. 
When the structure is not strictly rigid a purification pro- 
cess, based on the mean square method, makes it possible to 
"rigidify" it at the price of some approximations and a few 
more measurements. Lastly, eventual additional masses, that 
are not parts of the structure, can be taken into account. 



measurement of the inertial constants of a rigid or flexible 

STRUCTURE OF ARBITRARY SHAPE THROUGH A VIBRATION TEST 

D. Engrand, J. Cortial 


INTRODUCTION /214* 

A ground test is performed using usual equipment to deter- 
mine the inertial constants of the structure of an aircraft, 
missile, or of any other body whether flying or not. In many 
cases, the structure may be considered rigid, but potential 
deformations inevitably occur, and they are sometimes quite 
large. In R. Kappus's report [1] , the general case of a de- 
formable structure is first considered, and the resolution is 
obtained directly by methods of approximation (least squares) , 
by taking into account the spurious masses introduced by sus- 
pending the structure. 

As a user of this method, and we are especially addressing 
other potential users, we have considered the problem in a 
somewhat simpler form and we sometimes offer a different view- 
point than that of the author: in the first phase, the struc- 
ture is assumed to be perfectly rigid, suspended to a rigid 
frame using springs (linear) without mass. This somewhat 
simplistic scheme easily demonstrates that if we know the gen- 
eralized masses of the system (structure + suspension), and a 
convenient number of amplitude measurements, we can find all 
of the inertial characteristics of the structure: mass, co- 
ordinates at the center of gravity, moments and products of 

inertia in a given reference. 

In the second phase, the problem of a flexible structure 
is presented, by showing that it is possible, to a large ex- 
tent, to include the possible deformability of the structure 

if we know its first specific modes (for the type of boundary 
conditions imposed by the suspension) . 


♦Numbers in the margin indicate pagination in the original text. 
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In the last phase, corrections are made by accounting 
for the "spurious masses" introduced by the non-negligible 
masses of the suspension. This problem is presented very 
quickly in the appendix, for the simple purpose of providing 
a principle for eliminating the disturbances introduced by 
the suspension. Details of the operations required for this 
elimination are not indispensable for the understanding of 
this report; they were therefore limited to a few general 
characteristics. 

in '•losing this introduction, let us simply point out 
that we have we have attempted nothing more than to offer 
a summary, by sometimes adopting a slightly different view- 
point than that of R. Kappus, while retaining a system of 
notations very close to his. Specifically, the first three 
paragraphs are often quite different from the reference do- 
cument, while the last two are more direcly extracted from 
it. In the last paragraph (appendix) , we have not included 
the effective caclulation of spurious masses (matrix '0, 
so as to avoid needlessly cluttering a report whose purpose 
is simply to summarize the essential points of the reference 
document. 

I. - SIMPLIFIED DESCRIPTION OF THE PROBLEM FOR A RIGID STRUCTURE 
1.1 - Hypothesis 

-The center of gravity of the structure under study is 
assumed to be known. 

-The rigid structure S is suspended by a convenient number 
of negligible mass springs to an infinitely rigid and perfectly 
fixed frame. Let G„ be the position of the center of gravity 
G of S at rest. 

-The small motions of S are studied in the vicinity of 
the position of equilibrium, in a system of ortho-normed 
axes Goxyi 
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The solid thus suspended possess six degrees of freedom 
at the most. The motion will be described by a vector where 

x, y, z are the coordinates in a 0 x/7, 
and Ox, Oy, Oz are the rotations of S 
which are assumed to occur around cx, Gy, ci 
(hypothesis of small motions) . 1 is 

an arbitrary length intended to render 
the components of x. homogeneous. 

1.2. - Establishing The Equations 

The kinetic energy of the system 

in motion is expressed: 


G.M) 


xll 

% 

til 

h 

xll 

0 . 


2 r - m[/(C)|* + H.U1 


with: 



\Sr“ 

iU 

r<\ 1 

r 

no =. 

/ 

/ 


z 


LO 


and: 


hx ~ hr — ‘xz 
— hr hr — hi • 
.“hi —hi hi. 


( 1 . 2 . 1 ) 


( 1 . 2 . 2 ) 


(1.2.3) 


S . 



(tensor of inertia of S in G) 

The "linearized" expression (1.2.1) (i.e. omitting the 

terms whose order is greater than 2) is: 

2 T = m(x 2 + y* +■ *“) + t hr t hz^t ( 1 . 2 . 4 ) 

- 2 hrV\ - 2/ «V* “ 2 hAh- 
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If we let U‘be the quadratic form representing the 
"linearized" function of force from which the elastic sus- 
pension stresses are derived, the Lagrange equations are 
expressed: 

a, 2 . 8 ) 

m/ 2 O 0 0 0 0 

0 I 7Y O — I xr 0 — l Y2 

0 0 ml* 0 0 0 

0 /*). 0 I xx 0 — I xz 

0 0 0 0 m! 2 0 

0 ly 2 0 l X 2 0 I 22 
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(1.2.5) 


The system (1.2.5) may therefore be expressed in the 
conventional form: 

.of + ic/- o (1.2.6) 


with x being the stiffness matrix of the suspension system and /215 

IC the matrix of inertia explained in (1.2.5). « ■: i are the 

matrices defined as positive, symmetrical. We know that ma- 
trix .h-i .v. is, under these conditions, always diagonisable. 

Let Q be the matrix of its specific vectors. Matrices Q r . n 
and q t x<? in which 0 T is transposed by Q, are the diagonal 
matrices (orthogonality relationships of the specific forms) , 
and (1.2.6) may be expressed in the form: 

i/y + tfy Y «= o (1.2.7) 

with |* - 0 T .H0, Y = 0 t .K 0 , if = 0 -'.K, -, X 0 
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mu is the matrix of generalized masses, gamma the matrix of 
generalized stiffnesses, corresponding to the normalization 
selected for the specific vectors of • *.— 1 *• 

If we know Q, we see that it is possible to know , 
and therefore the inertias of S, by inverting according to 
relationship: 


•a. = (o’) Vo *. 


( 1 . 2 . 8 ) 


Generally speaking, a vibration test makes it possible 
to determine matrix mu. The operation for obtaining Q often 
presents more difficulties, as the structure is rarely rigid. 
Furthermore, the suspension springs are rarely without mass 
(or negligible mass with respect to that of S) , and the posi- 
tion of the center of gravity G of S is not always easy to 
determine on planes. 

We shall see in the next few paragraphs that it is not 
necessary to know point G a priori, and that it is possible, 
using a few approximations to overcome the experimental dif- 
ficulties. 

II - GENERAL EXPRESSION OF THE MATRIX JC 

(The position of the center of gravity G of S is not 
given) . 

1.1. Equations 

-Let 0 be any point of the structure. Its position 0„ 

C » 

at equilibrium will be adopted as the origin of a fixed 
reference oc*y* in which the motion of S will be described. 
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We will define a second reference oj-^. correlated with 
S, and which, at rest, is indistinguishable from ocxy*. 

The kinetic energy of the system is given by: 


■ - f \v (M )] 2 


dm (M) 


(II. 1.1) 


where M represents a running point or the solid, and we have: 


V (At) =-• V (0) + n A OM 

(ii being the vector for the rotation speed of S) 
We therefore have: 

(H.1.3) 2 [V (Atjj 2 dm « j [ / (0)] 2 dm + 

J s J 5 

+ 21/ (O).(U A ON) dm + 

r J : 

+ | ('ll A 0M).(ll A OAO d’". 


Noting that: 


J [/ (0)] 2 cfm = m[/(C)] 2 

! / (0). (Cl A OM) dm = 

Js . r 

= 2 (/ (0) A fl). I OM dm 


(II. 1.2) 


(II. 1.3) 


(II. 1.4) 


(II. 1.5) 


2m (V (0) A 0).OG 


f (H “ OAt). (ii A OM) 


rfm *= 12. J (0)11 


(II. 1.6) 


3 '.O) being the tensor of inertia in 0 of S, we may express: 
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2T = m [ / (O)] 3 + 2m (V (0) A ii) (II . 1 . 7 ) 

. O C + Cl. 3 (0)0. 

Note also that, is a vector of constant length 
involved in the me ion of the solid. We have: 


(n. 1 , 8 ) d"c = oo c + o c d 


oo, 


+ + « A O c O E 


(II. 1.8) 


with: 


g £ : • position of G at equilibrium (unknown) °* 

o„ • 

« - geometric rotation vector <•> ' 

If we explain the expression (II. 1.7) for kinetic energy 
by omitting the terms of order greater than 2, we derive: 


- m (x 2 + y* + /) + 2m(yf),—xil y )x a + 2 %V« 

+ 2m - *>, ,)y c + fy: . ~ 2 <Wxz ( 1 1 . 1 . 9 ) 

+ 2m (x« y — y l) x ) l a + *fyzz ~ 2,i * f Vxr 


With the following notations: 


x,y, z 

x G'y G ' z G 


Coordinates of 0 ( 

Coordinates of G (* n reference 0,-xyi. 

E 


Lagrange's equations may therefore be expressed: 


(II. 1.10) 

dtU<j,/ 


The first member having the following expression: 
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(II. 1.11) 


(11,1.11) </( \ 

j » m *’ + (''i./c — 
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pfl ' f : '.v.v' 1 , - Or'V - 0,0 + n ’ (*yo - y*c) 



?T\ .. , ’ 

i'i) ) ~ lyr >> ' ~ 0* 11 . + nl («, ~ * r c) 

dt ( 

d ( 

dt\ 

r n = - WV-W'x + m (y*o - *yc)- 

r'V 


As in the first case, the linear system forming small 
motions is expressed: 


Jt* 4- ax = o. 


(II. 1.12) 


With, in conformity with (II. 1.11): 
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(II. 1.13) 


(II. 1.14) 
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(Remember that 1 is the arbitrary length that we can 
choose according to the dimensions of the structure/ and 
that it does not serve to render the equations homogeneous) . 

Therefore, if we can return to matrix^- from matrix mu 
(generalized masses) , all of the inertial characteristics 
of the structure S are known: tensor of inertia in o and 
G coordinates in reference 0;^, and thus also the tensor of 
inertia in G in the same system of axes, and total mass m of S. 

III. - OBTAINING THE MATRIX FROM THE SPECIFIC FORMS Q 

The hypotheses are still those of the preceding paragraph. 
II. 1. - Preliminary Remarks 

-A point 0 of S being selected arbitrarily, the motion 
will be described by six generalized parameters that we will 
use -below as amplitudes of harmonic motions. These six para- 
meters form the components of the vector ; 


X 1 

z 1 
0* 

/ I 

L fj .J 


(III. 1.1) 


(vector a still relates to point )). 


Still under the hypothesis of small motions, any point 

P of S, of coordinates X , Y , Z , in reference oj-ri' will 

P P P _ 

have, in harmonic motions, a displacement given by: 
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We see that if we know x, 
rotation 
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we can determine the motion of 0. 


Therefore, in principle, if we can measure six independent 
amplitudes (two in each of the three directions Ox, o'y, 61 in 
the different points) it is possible to determine the motion 
of 0. 


in practice, we will do N measurements, for each of the 

six modes; at least two of the measurements will be performed 

in each of the directions d*. o'y, 61 Thus, for the k th mode, 

we will perform N measurements in N 1 points ,P N’ 

(N' is less than or equal to N) , with a maximum of three 

measurements (independent) made on each point along the three 

axes of reference o7~y7. These measurements will be stored on 

ad,,, vector: 

(k) 



(III. 1.3) 


we might of course have in column 5^, for a given point P^, 

the three components x p , y , z , or only two of these 

i *i 

instead of one. It should be stressed that the measuring 
points must be the same for the six modes. 


We may therefore express, in a similar manner to (III. 1.2); 

(III. 1.4) 
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i”"r, t 'o o 

*P. ° - x v> > Y n* O 0 


(matrix T having N lines and 
columns) . 

Note that matrix T does 
not depend on the measuring 
points selected i 


III. 2 - Obtaining the Specific /217 
Forms 


The structure being excited 
in the appropriate manner on 


its k mode, N measurement are performed which give the vector d (kr 

The specific form t will therefore be known if we are 

w 

able to solve the system: 

Th (k) - d (ky (III. 2.1) 

This system is generally superabundant. We perform its 
approached resolution using the least squares method. 

Resolution of the System (III. 2.1) 

If we call u. (j =1, ...., 6) the components of V ; and 
(i = 1, .... N) the components of d (ky the system may be 
expressed, with the convention of adding the silent indices: 


tv u, = v r 


The method consists of minimizing the quantity: 

1*1* - (fy "/- v ) (TV, 

by cancelling its partial derivatives with respect to 
u ^ ( = 1, ...., 6) 


(iix.2-2) 


(III. 2.4) 


3 




ORIGfr,'/! 

OE POO* QU* 


.$i v 


5!ff! 2 


C?u. 


f*U, 




or: 


"•(^“p •• ‘O = °- 


(III. 2. 5) 


If we vary t from 1 to 6, we obtain the new system: 


(III. 2. 6) 

r T Th m = T T d w . : 

Therefore: 

V,-( T'T)-' 

(II. 2. 7) 


(approached solution) . 

As the values of V> are given for the six modes by (III. 2. 7), 
the matrix Q is known and we can return to matrix .‘t. if W e know p, 

III. 3 - Determination of ^ 

This is obviously done experimentally, using in particular 
the "complex power method" [2] or the "displaced frequencies" 
method [3]. R. Kappus [1] recommends several measurements using 
several methods in order to be able to cross-check the results 
m order to be as accurate as possible, as the generalized mass- 
es are sometimes difficult to measure accurately. 

With the hypotheses of a perfectly undeformable solid, a suspen 
sjon without mass, and an infinitely rigid frame, we thus see 
that it is easy to return to the inertial characteristics of 
the structure, by starting with a somewhat modified ground 
test 13] . This case is clearly somewhat too ideal, and we 
are left with considering cases where the solid is deformable. 
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by taking into account the "spurious masses" and the possible 
deformation of the suspending frame, 

IV. - CASE OF THE DEFORMABLE SOLID 
IV. 1 - Hypotheses 

We are still assuming that the suspension is without mass, 
and that the frame is infinitely rigid, and perfectly fixed. 

In this paragraph, only the hypotheses made unti present of 
the undeformability of the structure S is discarded, and we 
are assuming that this structure is capable of underoing small 
elastic deformations. 


The gound test provides in this case a certain number v of 
specific modes (in f inifite theory) , wich will encompass both 
both the solid motions and structural deformations. Further- 
more, matrix .n can be determined only be means of the specific 
modes corresponding to the motions of pure solids. It is there- 
fore necessary to "rigidify" the structure, whenever possible, 
either by strengthening it mechanically during the ground test, 
or artifically, by looking for the specific modes of a fictive 
structure S' (attached to an identical suspension) offering 
the same geometery and the same mass distribution as the real 
structure at rest, but which should be undeformable. The meth- 
od proposed by Kappus [1] , and successfully applied, t rnsists 
of effecting a linear superposition of the specific modes re- 
corded, by assigning participation coefficients to them that 
are selected in such a manner that this superposition leads 
to a motion (fictive) that is as close as possible to the 
motion of a solid. 

In summary, this is a smoothing operation. The partici- 
pation coefficents are obtained using the least squares method. 
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IV. 2 - DESCRIPTION OF THE METHOD 


OK’ . 

OF POG . 


A number v(w v:.e> of specific modes are recorded during 
the test, each mode being represented by N measurements selec- 
ted as discussed in the preceding paragraph, and arranged in a 
column d (k) ' (k being the index of the mode). These columns, ar- 
ranged side-by-side, form a matrix V that will be the "measured 
modal matrix" (dependent upon the measuring points) : 

v lk = a - i th component of d (k) . (IV. 2.1) 


Note that V is a rectangualr matrix consisting of N lines 
and v columns. We still have to find a conversion enabling 
us to pass from the measured V matrix to the T matrix of the 
corresponding "rigidified" body. This is accomplished by 
looking for the participation coefficients -V; such as: 


\ - T „ 


(IV. 2. 2) 


(with the convention of silent indices) 
ving the system of matrices: 


This involves sol- 


r\ = T 


(IV. 2. 3) 


where A is a matrix consisting of v lines and six columns, 
which means that we have 6 <v equations for 6 v unknowns , 
and therefore a superabundant system. 


An approached resolution of this system can be easily 
obtained using the least squares method. As for (ill. 2. 6), 
this leads to the resolution of the system: 


V T V\ = V t T. 


(IV. 2.4) 
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As the columns of V are independent, as far as specific 

T 

modes are concerned/ V V is a constant matrix and we have: 

A ,(/Vr‘ V T r. (IV.2.5) 

While doing the calculations, it is of course useful to /218 

estimate the error introduced by the "least squares approxi- 
mation", although it is the best possible method (reference 
II] , pp. 12 and 13) . 

We then have to return from matrix p (that now has v ; 
lines and v columns) to matrix .*• 

IV. 3 - CALCULATION OF MATRIX 
a) Correlation of Orthogonality; Matrix M 

In any specific motion of the system, of angular velocity 
w, represented by a vector d of which N components are the 
N amplitudes measured, the vibrations of the system are govern- 
ed by a system with the expression; 

^Md-Kd (IV. 3.1) 


M and K are the squares matrices of order N, representing 
the mass distribution and rigidities of the system (S + sus- 
pension) . In theory, these matrices may be explained on the 
basis of the total energy of the system in any vibratory mo- 
tion. 


In fact, they will be precious intermediaries, but it will 
not be necessary to explain them. 

-The correlations of orthogonality of the specific vec- 
tors will be expressed in this representation; 
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(IV. 3.2) 


V t MV • a 


1 1*11 w 
I 1*22-: 

I 0 V 


For a give structure M, M depends only upon the selected 
measuring points. 


b) Correlation Between M and 


For any virtual harmonic motion, of angular velocity w 
in conformity with the undeformability of the "rigidif ied" 
structure, defined by: 


a * 


v/f 

°r 

z'n 

o; 

y'n 


(IV. 3. 3) 


the kinetic energy E* is expressed: 


s' = 



(IV. 3. 4) 


The corresponding virtual "measurement vector" is 
expressed: 


Tn 


(IV. 3. 5) 


and the kinetic energy may also be expressed: 


or 


(IV. 3.6) 

S' = V* <)* r AW* 

2 

s' - \,>’*i’ T r T A :ri*. 

r\ 

(IV. 3. 7) 

vuW. ~~ 
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Through identification between (IV. 3. 4) and (IV. 3. 70, we 
therefore obtain: 


n T t MT 


(IV. 3 .8) 


By correlating (IV. 2.3) and (IV. 3.*), we therefore have: 

.* - \ T V T AfKA. (IV. 3. 9) 

(We should still not forget that this equality is not 
perfectly accurate due to the approximation made in the cal- 
culation of A).. 

By using the correlation of orthogonality (.IV. 3. 2), we 
finally obtain: 



(IV. 3.10) 


Note: This method of resolution is valid only if the 
deformations of S occur only in the form of spurious motions, 
which is often the case. It constitutes the main point of 
the method developed by R. Kappus II] , and is also its main 
advantage. 

APPENDIX - ELIMINATION OF SPURIOUS MASSES 


Al. - 


For the report to be complete, it is necessary to bring 
up this problem, because it is not always of the secondary 
importance that we have described in the preceding paragraphs. 

The expression" spurious masses" designates here all 
masses contributing to the motion and that do not belong to 
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the structure S itself. We are dealing with suspension rods 
and, possibly, the frame when the latter is not rigid enough 
(situation that is good to avoid) . Strictly speaking, the 
inclusion of these masses means that additional modes will 
have to be taken into consideration, but we will assume that 
these be only of a very local nature and that they have spe- 
cific frequencies that are high with respect to those measured 
during the ground tests. 

A2 - DESCRIPTION OP THE METHOD [1] pp. 6 to 8) 
a) Measurements 

N measuring points should be provided and distributed 
over the structure and suspension. N q points (N q is less 
than or equal to N) will be selected on the structure, N - N Q 
on the suspension (possibly on the frame also) . We will re- 
cord •< specific modes (S is always assumed to be deformable, 
and therefore •' ,6) - Therefore: 


The modal matrix V is relative to N measurements, as well 
as matrices M and K. We still have the correlation._r)f.„-or±ha^. 
gonality : 


According to the assumptions made, f* is disturbed by 
masses foreign to S. The purpose of the following is to de- 
fine a matrix rnu Q that correlates only to S, using a process 
of elimination of the spurious masses. 

b) Resolution of the Problem 

It consists of subtracting from M a matrix M' representing 
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the effects per unit of mass of the suspension rods (or 
possibly the frame) , and of looking for a modal matrix V q 
representing the measurements that we could have performed 
if the suspension was without mass. This matrix V q verifies 
a theoretical- relationship expressed: 

v a = (a2 

(Omega o being the unknown matrix (diagonal) of the specific 
angular velocities of the system without spurious masses). 

The method proposed by R. Kappus consists of setting: 

V = Vr (A2 

o 

R being an unknown least squares matrix of order v , 
(A2.4) becomes: 


(M — M‘) vr = kvrq 0 2 . 


(A2 


This is a superabundant system that we solve again 
using the leas t sq uares method: 

\v T mv- v t m v)Rm (v T kv) mV 2 - 

Accounting for the correlations of orthogonality, and 
by setting v' T M/ = r or assumed to be known), we obtain: 

(AS 

(remember that omega is the matrix of the measured specific 
angular velocities). 


The system (A2.6) has a conventional equation: 




jr = ir fib -2 

where I and J are the symmetrical matrices defines as posi- 
tive. We know that thesolution R of such a system is an 
orthogonal matrix in terms of J. In other words, the matrix: 


a a = R T JR - R 


is a diagonal matrix. As R is thus determined, we have V Q = 

VR, and we now only have to truncate the matrix V Q by removing 

its lines corresponding to the measuring points foreign to S 

itself. Let V be this truncated matrix. It is also useful 
OS T 

to define a matrix M q (N X N) such that F- Q = V Qg M 0 V 0S • 

We then simply have to finish solving the problem by 
replacing N by N , V by V Qg , M by M q , R by f* Q in paragraph 
IV. 

c) Obtaining Matrix M' 

This matrix (symmetrical) must be determined by calcula- 
ting on planes ([1] pp. 14 to 17), which may be relatively 
simple. 

CONCLUSION 

We have seen that, even in the most unfavorable cases, 
it is possible to determine the inertia constants (tensor of 
inertia, mass, center of gravity) of a structure by doing a 
few approximations on a simple ground test. Generally speak- 
ing the procedure to follow is: 

1. Eliminate the spurious masses, if necessary. 

2. Solve the problem according to the method in paragraph 
IV if the structure is deformable, or according to that of 
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paragraph III if it is found to be undeformable . 

In any case, the resolution is always carried out using 
an exisitng automatic computer program at the ONERA computer 
center. Numerous applications [3] have made it possible to 
evaluate the method with considerable accuracy. 


This document was submitted on June 4, 1970. 
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